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We investigate the morphology of diblock copolymers in the vicinity of flat, chemically patterned 
surfaces. Using a Ginzburg-Landau free energy, spatial variations of the order parameter are given 
in terms of a general two-dimensional surface pattern above the order-disorder transition. The 
propagation of several surface patterns into the bulk is investigated. The oscillation period and 
decay length of the surface g-modes are calculated in terms of system parameters. We observe 
lateral order parallel to the surface as a result of order perpendicular to the surface. Surfaces which 
has a finite size chemical pattern (e.g., a stripe of finite width) induces lamellar ordering extending 
into the bulk. Close to the surface pattern the lamellae are strongly perturbed adjusting to the 
pattern. 



I. INTRODUCTION 

The bulk properties of diblock copolymers (BCP) are 
now well understood These long linear macro- 

molecules composed of two incompatible sub-chains, or 
blocks, cannot phase separate because of the covalent 
bond between them. This connectivity, together with the 
incompatibility between the two blocks, gives rise to the 
appearance of microphase separated phases. The state 
of segregation is controlled by Nx and /, where x is the 
Flory parameter, TV — Na + Nb is the total number of 
constituents monomers per chain and / is the fraction of 
the A block, f = Na/ {Na + Nb). For large enough Nx 
one of the ordered phases, such as the lamellar, hexago- 
nal or cubic phases is preferred, depending on the degree 
of asymmetry /. 

Less understood is the interfacial behavior of copoly- 
mer melts near solid surfaces or at the free surface with 
air. Surface phenomena of BCP may enable creating and 
controlling technologically important devices of charac- 
teristic size comparable to the wavelength of light. As 
examples we mention waveguides, light-emitting diodes 
and other optoelectronic device, anti-reflection coating 
for optical surfaces ^ and dielectric mirrors 0. 

The presence of a wall in a BCP system leads to new 
energy and length scales, depending on the specific chem- 
ical interaction of the polymers with the surface. In a 
semi-infinite system in contact with a single planar wall, 
the morphology near the surface can be very different 
from the bulk morphology. Fredrickson [^j has consid- 
ered BCP in contact with a surface having a uniform 
preferential adsorption to one of the two blocks. Above 
the order-disorder transition (ODT), where x < Xc (Xc 
is the critical point value of x above which an ordered 
phase appears), he used mean-field theory and found 
that the order parameter (being the concentration differ- 



ence between the two blocks) has decaying oscillations. 
He showed that the oscillation periodicity depends on x, 
and tends to the bulk lamellar periodicity as x ^ Xc- 
In the same X ^ Xc limit, the correlation length ^ of 
the oscillations was found to diverge. Further investiga- 
tions H] showed that the inclusion of higher order, non- 
linear corrections to the mean-field theory results in a 
non-diverging i^. For the same system cooled below the 
ODT, modulated sinusoidal behavior was found. In a re- 
lated work a Ginzburg-Landau free energy was used 
to describe the propagation of a surface-induced lamel- 
lar ordering into a bulk hexagonal phase. In the strong- 
segregation limit a lamellar region of finite thickness close 
to the surface becomes stable, provided that the surface 
field is larger than some critical value. 

The situation is even more complicated in thin films, 
where the distance between the two boundaries, associ- 
ated with the film thickness, is comparable to the period- 
icity of modulations in the bulk, and the surface induced 
morphology can be of different symmetry than that of 
the bulk. For a system taken in one of its ordered phases 
(below the ODT), the free energy has a local minimum 
when the spacing between the surfaces is an integer mul- 
tiple of the bulk repeat period. The mean-field behavior 
of BCP close to surfaces and for BCP films was calcu- 
lated ||ri[] using a method applicable in both the strong 
and weak segregation limits. It was found that confine- 
ment of lamellar phase BCP may lead to parallel layering, 
or in some cases even to a perpendicular arrangement of 
the lamellae. Self-consistent field theory (SCF) was used 
Jl2| , p^ to study the stability of these parallel, perpendic- 
ular and mixed lamellar phases in thin films of BCP. The 
latter phase consists of parallel lamellae near one surface 
and perpendicular lamellae near the opposite surface, but 
it was found to be unstable for symmetric A-B (/ = 1/2) 
diblock copolymers. 

So far, we mentioned situations where the surfaces have 
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a uniform preference to one of the two blocks. More 
complex, chemically patterned surfaces break the lateral 
translation symmetry. Different surface regions will now 
have a different preference for the A/B blocks, thereby in- 
ducing a lateral structured morphology near the surface. 
Very few works took into account this possibility of a non- 
uniform surface. In particular, Petera and Muthukumar 
have investigated the effect of a one dimensional 
sinusoidal surface pattern on BCP morphologies close to 
the surface in the weak-segregation limit, both below and 
above the ODT. 

In this paper we consider a BCP melt above the ODT 
near a surface, whose pattern is truly arbitrary in two 
dimensions, generalizing the results of Refs. A 
Ginzburg-Landau free energy is expressed in term of the 
polymer concentration is presented in Sec. In Sec. 
m we consider a melt close to one surface or confined 
between two surfaces whose chemical pattern has one di- 
mensional symmetry. Minimization of the free energy 
expansion gives rise to an Euler-Lagrange equation for 
the order parameter. A natural generalization to two- 
dimensional surface patterns is then considered in Sec. 
IV . We are able to give a complete description of the or- 
der parameter in terms of all the 5-modes of the surface 
pattern. Finally, conclusions and some future prospects 
are presented in Sec. 



II. THE MODEL 

The copolymer melt is described by the order param- 
eter (j){r), defined as 0(r) — (/)yi(r) — /, the difference in 
local A monomer concentration from its average value. 
Hereafter we restrict the treatment to the symmetric 
/ = 1/2 case, following the same coarse-grained free en- 
ergy as was used by Fredrickson and Binder p6| , pT| : 
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Where ks is the Boltzmann constant and T is the tem- 



perature. The other parameters are: 




go « 1.9456/ 


(2) 


T = 2pN (xc - X) 


(3) 


Xc ^ 10.495/A^ 


(4) 


h = \.bpc^{Rl)/ql 


(5) 



The fundamental wavelength of the system, qoi is ex- 
pressed by the radius of gyration of the chains. The 
chain density p is equal to for an incompressible 

melt, and u/ p and c are of order unity. More details can 
be found in Ref. ||l^ and extensions for asymmetric BCP, 
/ 7^ 1/2 are possible as well. The use of (|lj) limits our 
treatment to a region of the phase diagram close enough 
to the critical point where the expansion in powers of (f> 



and its derivatives is valid, but not too close to it, be- 
cause then critical fluctuation effects may be important 

This and similar types of free energy has been used 
to describe bulk and surface phenomena in amphiphilic 
systems , diblock copolymers p|,^, p^,p^j2l| , Langmuir 
films and magnetic (garnet) films |23|[~^The 0^ and 4>'^ 
terms appear in the usual Landau expansion. The added 
(/)V^(/) and 0V^V^(/) terms compete to produce modulated 
phases below the order-disorder temperature. This free 
energy describes a system in the disordered phase (0 = 0, 
/ = 1/2) for X < Xc, and in the lamellar phase for x > Xc- 
The q = qo mode goes critical first, and the lamellar 
phase is described by (/> = 0^ cos(qo • r), of repeat period 
do = 27r/go. This single-mode approximation is accurate 
to order {x~ XcY^^ and can be justified near the critical 
point Far from the critical point higher harmonics 
are needed to describe the lamellar phase. As the asym- 
metry in composition is increased, other ordered phases 
of hexagonal and cubic symmetries become more stable 
than the lamellar phase. 

As stated above, block copolymers exhibit complex 
surface behavior characterized by the strength and range 
of the interaction between the polymer chains and the 
surface, the typical size of chemical heterogeneities of the 
surface, and the distance between the two surfaces, in 
case of a thin film. 

The presence of chemically interacting confining walls 
is modeled by an added short-range surface coupling term 
in the free energy. 
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The vector Tg define the position of the confining sur- 
faces. The acj) term expresses the preferential interaction 
of the surface with the A and B blocks. For example, if 
a > then the B block (0 < 0)is attracted to the sur- 
face more than the A block (0 > 0). Control over the 
specificity of this surface term can be achieved by coating 
the substrate with carefully prepared random copolymers 
p^ , p5[ . The coefficient of the 0^ term in (ph , , is a sur- 
face correction to the Flory parameter x P,R|2^ . > 
corresponds to a suppression of surface segregation of the 
A and B monomers. 

We first consider systems in which the polymer melt 
is confined by a flat, rigid wall at y = 0, with the a;-axis 
chosen in the plane of the wall, and is translational in- 
variant along the z-direction. Extension to the system 
of two parallel surfaces located at y = ±L is straightfor- 
ward and will be considered later. The order parameter 
(j) is expected to vanish in regions where the interfacial 
interactions can be neglected. 
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recovering the value = of the bulk phase far from 
the surface. In the next section we find profile solutions 
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4>{x, y) for a BCP system at temperatures above the bulk 
ODT. 



III. ONE DIMENSIONAL SURFACE PATTERNS 

For high enough temperatures, or equivalently, for 
X < Xc, the phase of lowest free energy is the homo- 
geneous disordered phase, with 0(r) = in the bulk. 
The presence of a patterned surface induces ordering in 
the copolymer melt. If the chemical surface composi- 
tion is uniform, one of the monomers, A or B, will be 
attracted to the surface, resulting in a parallel orienta- 
tion of the lamellae (a perpendicular orientation of the 
chains). If the pattern is modulated, say sinusoidally, 
then different blocks are attracted to different regions of 
the surface. We will start with this case of a semi-infinite 
system bounded by one rigid, flat surface, and then pro- 
ceed to describe thin film systems between two surfaces. 



A. One patterned surface 

Consider the semi-infinite BCP melt at y > bounded 
by a flat surface given by y = 0. We assume a one di- 
mensional periodic surface pattern and write it in terms 
of the Fourier components of the surface field a{x) 



a{x) 
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where aq set the amplitude of the respective g-modes. 
The order parameter (/)(a;, y) satisfies the boundary con- 
ditions on the surface and approaches the bulk solution 
far from the surface. It is convenient to decompose ((> in 
terms of its g-modes in the a;-direction 



0(^,2;) =^/,(y)e^' 
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The requirement (0) leads to the bulk boundary condi- 
tion 



hm fqiy) = 
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The form (^) is substituted in (|i|). Above the order- 
disorder transition (ODT) temperature, the theory is sta- 
ble to second order in (p, and therefore the (j)'^ term is ne- 
glected. Using the explicit x-dependence of cj) in (^) we 
perform the x and z integration, yielding the free energy 
F: 

^ = / E { + '^'^o) fj; + hql [h - q^h) * + c.c- 



+ lh[fq[fq" -2q^f'; +q^fqy 

+ EK/9(0)+^«/9(0)/;(0))+c.r 



dy 
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where (...)* indicates complex conjugation ( c.c. ) oper- 
ation. A standard minimization technique is carried on 
and yields the governing linear ordinary differential equa- 
tion for the fimctions {fq} for y > 0: 

(r/h + (g2 _ q^)'^ + 2{ql - q^)/; + /;' - (12) 

This equation possesses four independent solutions in the 
form of an exponential e~''«^, with kq found from the 
characteristic equation 

{r/h + {q' - qlY) + 2{ql ~ q^)kl + fc^ = (13) 

Thus, with the semi-infinite geometry, the solution is 

fq{y) = Aqe-'^-y + Bqe-Ky (14) 



and 
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Each g-mode solution fq is characterized by two complex 
amphtudes {Aq.Bq}. From the solutions of Eq. ( |l^ ) 
wavevectors {kq} with negative real value, Re(fcq) < 0, 
are discarded, to comply with the boundary condition 
(pO|). Note that Re(fcg) is increasing monotonously as a 
function of q. A large value of Re(fcg) means short de- 
cay length, and hence the smallest surface g-mode decays 
the least. This behavior is demonstrated on Fig. 1 show- 
ing the q and x dependence of the real and imaginary 
parts of the wave- vector kq. For a fixed value of x ths 
decay length, proportional to 1/Re(fcg), decreases as q 
increases, while the wavelength 27r/Im(/cq) of the modu- 
lations in fq{y) ^ e^'^iV increases. 

The boundary conditions for the functions fq can be 
determined by considering the Euler-Lagrange equation 
for {fq} in the range that includes y = 0. In this case 
a term proportional to the Dirac delta function 5{y) ap- 
pears in (|l2|). There are two conditions relating fq and 
its derivatives at y = 0: 

fl{Q) + 2{ql-q^)fq{Q)^0 (16) 

^ + + 2(9o' - 9')/;(0) + /;(0) = (17) 

Recently, it has been found that surface states ex- 
ist even in the absence of a surface field a. This effect 
can be attributed to a loss of entropy close to the sur- 
faces. However, in our linear theory this does not hap- 
pen, and the copolymer response is proportional to the 
surface field a. The case where ctq is a non-zero constant 
and Cg^o — 0, corresponds to the special case of uniform 
interfacial interactions. The system exhibits a decaying 
lamellar layering of the polymers, with the B-polymer 
adsorbed to the surface if (Tq > 0. 

Close to the ODT, the complex wavevector fco can be 
approximated by 
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This expression shows that /o ~ e '''"^ ~ e has a di- 

verging characteristic length ^ oc (xc — X) ^ i while the 
oscillatory part has a wavelength slightly longer than that 
of the bulk lamellar phase, in agreement with the results 
of Fredrickson for chemically uniform surfaces The 
correlation length ^ diverges for small composition os- 
cillations because a linear theory is employed; addition 
of the (/)^ term in F would give a finite value of ^. As 
a result of the assumed short-range surface interactions, 
the periodicity and decay length of fq depend only on 
properties of the bulk, and not on surface details. 

Using the notation kq = + ik^ the real part of the 
form (|lj) can be rewritten as: 



2Re {fq) = {Aq + B; 



2\Aq + Bq\e-''iy cos{k'^y + a,) 



(19) 



where aq is the phase of the q-mode modulation. It de- 
termines the value of the fq solution at the boundary, 

The phase ao for the q — mode is found to satisfy 
the following relation 



tan ao 



(20) 



and therefore is determined by the degree of segregation 
X, but not by the pattern amplitude aq. A plot of aq as 
a function of q for several values of x is shown in Fig. 2 
(a). Far from the ODT point and deep into the disor- 
dered phase, X ^ Xc, we find that all g-modes have a 
phase angle aq = 0. As the ODT is approached, the 
q = qo mode retains its value but larger q-modes have 
a negative phase while smaller g-modes have a positive 
phase. At X = Xc this becomes a step function, with 
aq = 7r/2 for q < qo and aq = — 7r/2 for q > qo- The 
amplitude behavior is shown in Fig. 2 (b), for the same 
series of segregation values x- As the ODT is approached, 
the q = qo mode becomes critical first, with a diverging 
amplitude. 

An interesting limit occurs when ctq = 0, that is, the 
average surface interaction is zero (no net adsorption). 
No lamellar ordering parallel to the surface is expected. 
Indeed, the resulting checkerboard behavior is illustrated 
for a surface pattern chosen for simplicity to contain only 
one mode: (t{x) — aqCos{qx). Fig. 3 depicts alternating 
A-rich (white) and B-rich (black) regions. In (a) the de- 
cay length ^ is smaller than in (b) , because in the former 
case the surface periodicity is twice as large. The oscilla- 
tory behavior, characterized by lm{kq), has a very long 
wavelength, diverging as (Xc - x)"^/^ close to the ODT 
point. 

Usually, if no special measures are taken [ p4p5| ] , there 
is a net preference to one of the monomers: ao 0. 



The BCP morphology where the surface interactions 
were chosen to have both a non-zero average preference 
and undulatory character, namely a = ao + aq cos{qx), 
is shown in Fig. 4. A smooth crossover from surface- 
induced ordering at small distance to the bulk disorder 
occurs. The parallel lamellae resulting from the ao term 
persist farther from the surface than the bulges resulting 
from the aq term, as /q decays slower than fq. For a 
given ao, having a higher q-mode or reducing the modu- 
lation strength aq will enhance the lamellar features far 
from the surface. 



B. Two patterned surfaces 

Until now the BCP melt was assumed to be bounded 
by one surface at y = 0. In this section we extend our 
analysis to a thin-film system confined between two par- 
allel surfaces located at y = L and y = —L, shown in 
Fig. 5. When the distance 2L between the surfaces is 
comparable to the natural bulk periodicity, the two sur- 
faces interact via the BCP and the resulting film mor- 
phology ca n be v ery different from that of the one-surface 
case (Sec. [II A). However, the mathematical analysis is 
almost the same; one only has to apply different bound- 
ary conditions on the BCP order parameter (p. 

The surfaces at y = ±L are assumed to carry differ- 
ent surface fields of the form a^{x) — J2q cy^Q^'^^ ■ Only 
small modifications must be included to adjust the re- 
sults of the previous section. The same ansatz (|^) for (50 
is used. The functions /„ that minimize the appropriate 
x-averaged free energy d|) obey 



fq{y) = Aqe-'^-y + Bqe-'^'^y + Cqe'^^y + DqeKv (21) 



with {kq} given by the same relation (15). However, 
unlike the semi-infinite bulk one-surface case (|4|), both 
signs of the k vectors are used because the system is finite 
in the y-direction. In addition, repeating the procedure 
outlined for the one-surface case, gives the boundary con- 
ditions for fq-. 



f';{±L) + 2 {ql - q^) fq{±L) = 



2a^ 



4t, 



h h 

T2 {ql-q')f'^{±L)Tf'q"{±L) = 



(22) 
(23) 



We consider now several specific surfaces. In the first 
= 1 is a constant, while a~ = cos{qx) is purely si- 
nusoidal and average to zero, as depicted in Fig. 6. As 
expected, the B polymer (in black) is attracted to the 
upper surface, while the bottom surface exhibits modu- 
lated adsorption pattern. Although lamellar features are 
seen near the top surface, the overall apparent phase in 
the sample cannot be classified as such. The correspond- 
ing plots of the functions fo{x) and fq{x) are shown in 
Fig. 7 (same parameters as in Fig. 6). In general fq is 
nonzero even at the y — L surface, although the surface 
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does not induce modulations by itself, cr"*" = const. Thus 
modulations propagate from one surface to the other by 
the copolymer melt. This is an interesting observation 
which may have relevance in applications. It relies on 
the relative small thickness of the BCP film. 

The situation as depicted in Fig. 6 represents a com- 
petition between two mechanisms. The modulated pat- 
tern at the bottom surface induces a laterally modulated 
pattern of the BCP, while the top surface uniform inter- 
action induces a lamellar-like layering of the copolymers. 
As the modulated adsorption pattern strongly depends 
on the modulation wavenumber, so does the resulting 
morphology. This effect is shown explicitly in Fig. 8, 
where the top surface is uniform and the bottom is mod- 
ulated, for a series of q/qo values. The transition from 
a locally perpendicular (bottom patterned surface) to a 
locally parallel orientation (at the top uniform surface) is 
seen in (a), similar to the so-called T-junctions between 
grains of different orientations ||2l|,^. Similar behavior 
was found by the SCF calculation in Ref. ||l5[] . 

Figure. 9 (a) shows the spatial dependence of the BCP 
order parameter when the two surfaces contain only one 
g-mode and are patterned in phase with each other (sym- 
metric arrangement), but with opposite signs, cr^ = 
±aqCOs(qx). The copolymer patterns create a perfect 
checkerboard arrangement and are related to each other 
at the surfaces by an interchange of monomers A-s-^B. 
The surface pattern (||) contains only cos(gx) terms. A 
generalization that includes sin((7x) sinusoidally varying 
modes is straightforward. In this case the patterns at 
the surfaces can be out-of-phase with each other. Fig- 
ure. 9 (b) shows such a morphology, for = Uq cos{qx), 
— CT^sin(ga;), where there is a 7r/2 phase shift be- 
tween the two surface fields. The perfect checkerboard 
arrangement of 9 (a) is now distorted to accommodate 
this phase shift. 



IV. TWO-DIMENSIONAL SURFACE PATTERNS 

So far we considered a melt in contact with a surface or 
confined between two surfaces of one dimensional symme- 
try. In our approximation, like in any linear response the- 
ory, there is no g-mode coupling proportional to ■ 
This fact allows us to go further and introduce a two- 
dimensional generalization of the surface pattern, which 
so far was taken to be independent on z. The surface now 
is assumed to carry a chemical pattern (t{x, z) which can 
be written as: 

a(x,2)= ^ a,,,,,e'(«--+«^^) (24) 
qx,qz 

The "linear response" function is then 

50(x,y,z)= ^ fq^^qAyy^'-'''^'^'^ (25) 



Because / and a are real functions, f-q^^-q^ = f*^ 
and similarly for a. Using the above form it is possible 
to carry out the integration of the free energy in the x — z 
plane. Denoting {■■.)xz as the average in the x — z plane, 
it can be checked, for example, that 

{cl>V^)x. = fq^,q^ {f'qlq^ - {ql + ql) fq^,qS (26) 

Defining = {qx,qz) and performing the free energy 
minimization with respect to f*^ , the functions fq^^ — 
fq^^q^ obey the same master equation ([l^ ) that fq pre- 
viously obeyed, with the only change that q^ is replaced 
by q^. For a BCP in contact with a single surface, the 
appropriate boundary conditions are: 

f'q[i0) + 2(^ql-q^^fq^^{0)^0 

^ + ^fi« (0) + +2 {<lo ql) /;, (0) + /;'(0) = (27) 
The solution for fq^^ is analogous to (p^, 

/,„(2/) = ^,„e"S^ + B,,,e-''V (28) 

kl=q^-l+^{Xc-xf' (29) 

Having found the response of the polymers to the sur- 
face modes aq^^ , one is able to deduce the concentration 
profiles for any given two-dimensional surface pattern. In 
order to illustrate this, we take a system of chemical affin- 
ity in the shape of the letters "BCP" on the y = surface 
[Fig. 10 (a)], and calculate the polymer concentration in 
the planes parallel and above it. All sizes are expressed 
in terms of do, the lamellar fundamental periodicity. The 
shape of the letters continuously deforms as one moves 
away from the surface. Contour plots corresponding to 
planes parallel to the x — z surface and separated by a dis- 
tance {n+ \/2)da, for integer n, are approximately given 
by an A ^ _B interchange of monomers. This is the char- 
acteristic distance at which the polymers flip. Note that 
Fig. 10 (c) and (e) are approximately the inverse image 
("negative") of (b), (d) and (f). The original features are 
completely washed away as the distance y from the sur- 
face is further increased. In our case for Sdo ^ y ^ 6do 
where the surface pattern size was roughly do. 

Figure 10 also illustrates circumstances where a cer- 
tain surface pattern is transferred via the bulk BCP to 
another, distant surface. It may be important to know, 
for example, if the contrast of the distant image can be 
experimentally detected. This reduction of the contrast 
is clearly seen by comparing Fig. 10 (b) and (f), and can 
easily be calculated from our expressions. The lamel- 
lar order created parallel to the edges of the letters in 
Fig. 10 (b) is the result of the undulatory nature of the 
block copolymers: order extending perpendicular to the 
surface induces order in the direction parallel to it. 
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The copolymer melt can follow the surface pattern 
when its size is larger than the polymer length-scale do. 
The effect of reducing the size of the surface structure is 
seen as a blurred morphology in Fig. 11 (a), where the 
"BCP" pattern was chosen to have dimensions Ado x 4(ioi 
compare to 20c?o x 2Qdo of Fig. 10 (b). The effect of rais- 
ing the temperature (further away from the ODT) is seen 
in Fig. 11 (b). It is similar to Fig. 10 (b), only that the 
temperature is higher, Nx = 9, and the lamellar features 
along the edges of the letter are less prominent. Again, 
using our order parameter expressions one can quantify 
the g-mode spectrum and contrast of the distant image, 
as a function of the original surface pattern a{x, z), tem- 
perature and distance from the y = surface. 

In a thin film, creation of truly three dimensional, 
complex morphologies between the two surfaces can be 
achieved by using only one-dimensional surface patterns. 
As an example we choose a simple sinusoidal pattern 
on each of the y = ±L surfaces, rotated 90 degrees 
with respect to one another: a~ = cos{qx + qz) and 



r+ - 



cos{qx 



qz). In Fig. 12(b) the resulting morphol- 
ogy in the y = mid-plane is shown and is a superposi- 
tion of the two surface patterns. Because y = is a sym- 
metric plane, the pattern has a square symmetry. More 
complex patterns can be created at different y planes. 
In Sec. Ill we showed that the q = mode of the 



surface pattern is the slowest decaying mode, resulting 
in a lamellar layering parallel to the surface as y ^ cxd, 
no matter what the surface pattern is. We demonstrate 
this in Fig. 13, where in (a) we choose a simple one- 
dimensional structure in the shape of a stripe of width 
do. Inside the stripe of width do, a{x,z) = 0.5 while 
outside it, the surface area is neutral, a — 0. Thus, the 
B-polymer is preferentially adsorbed onto the stripe. The 
order parameter contour plot in the x — y plane is shown 
in (b). It can be seen that the "surface disturbance" is 
enclosed with alternating lamellae. The distorted lamel- 
lae close to the y = surface appear curved, and slowly 
fade away as the distance from the surface is increased. 

A different scenario is presented in Fig. 14, where in- 
side the stripe of thickness do: = 0.5 as above, but 
outside the stripe the surface is not neutral: a = —0.5. 
We find that the adsorption on the surface is quite differ- 
ent than in Fig. 13. Far from the stripe, the A-polymer 
is adsorbed onto the surface and induces stacking of the 
BCP in a direction parallel to the surface. Close to the 
surface perturbation (the stripe) the behavior is altered 
as the lamellae are strongly deformed in order to optimize 
their local interaction with the surface stripe. 



V. CONCLUSIONS 



We have employed a simple Ginzburg-Landau expan- 
sion of the BCP free energy to study analytically the 
confinement effects of block copolymers between two pat- 
terned surfaces as well as the interfacial behavior of a 



BCP close to a patterned surface. Our approach consists 
of finding the governing equation for a presumably small 
perturbation to the bulk order parameter, by retaining 
second-order terms in the free energy. This approach can 
be justified in the vicinity of the critical point. Above the 
ODT it gives rise to a simple linear equation with fixed 
coefficients P,p^. A generalization to two-dimensional 
surface patterns is presented, where a complete spatial 
description of the polymer concentration is given in terms 
of an arbitrary surface pattern. However, this approach 
applies to systems below the ODT as well, where a lin- 
earization is to be taken around an ordered phase |^,|lj . 

The assumption that the surface interactions are 
strictly local means that the length scales of the polymer 
morphology are determined by bulk properties. More- 
over, each of the surface g-modes in a{x) = aq cosqx 
gives rise to a corresponding mode fq cos qx in the local 
polymer concentration (f){x,y). This "response" mode is 
characterized by a single wavevector kq . The wavevector 
kq is determined by x and the surface wavenumber q. In 
Fig. 1 we show the dependence of kq on these parameters. 
The high g-modes of the surface pattern a{x, z) decay 
more rapidly than those of low q. For high g-modes of 
characteristic length scale much smaller than the poly- 
mer chains do, the BCP melt cannot follow the surface 
modulations, and feels just the average of those modula- 
tions (which is zero for q > 0). This dependence of kq on 
q and x is very similar to the results found by Petera and 
Muthukumar [ p^ using a different free-energy functional. 

Moreover, we generalized surface patterns to any two- 
dimensional patterns as can be seen in Fig. 10. Even 
within a mode decoupled (linear response) theory, many 
interesting effects follow for a single surface as well as for 
films confined between two surfaces. Tuning a few sur- 
face parameters can lead to controlled micro-structures 
of the BCP film. 

For a BCP melt in contact with a homogeneous sur- 
face, a decaying lamellar order appears. The phase ao 
of these sinusoidally damped oscillations obeys tan ao = 
qg/ y^r/h, and it is independent of the surface pattern 
amplitude ctq |]i7| . For high temperatures all q- modes 
have the same phase = 0. As x Xc, the q = qo 
retains this value, while higher g-modes tend to — 7r/2, 
and lower q tend to 7r/2. At the same limit the q — qo 
mode gets critical first, with a diverging amplitude. 

Our expressions for the spatial dependence of the order 
parameter on a general patterned surface gives a com- 
plete description of the system, and allows for the calcu- 
lation of free energy, pressure, etc. It may also help in 
tuning the required distance between the two surfaces in 
various applications. Using a strong enough surface field 
or fixing the conditions close to the ODT, one can hope to 
transform a pattern from one surface to the other surface. 
We also demonstrate in Fig. 12 how the superposition of 
simple one dimensional patterns can bring about a three 
dimensional behavior in a thin film system. A desired 
complex phase can then be achieved by tuning the Flory 
parameter and the relevant distances. 
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Possible extension to this work will be to calculate the 
phase diagram of the L± phase (a confined lamellar phase 
where lamellae are perpendicular to the confining sur- 
faces) vs. the L|| phase (where the lamellae are parallel 
to the surfaces), by calculating the surface contribution 
to the bulk free energy in F (^. In the weak segrega- 
tion limit this contribution is important and may lead to 
a completely different diagram than that of the strong 
segregation regime. 
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• Fig. 1: The real (a) and imaginary (b) parts of 
the wavevector kq as a function of the modulation 
g-mode and the Flory parameter Nx- For values 
of Nx close to its critical value, Nxc = 10.495, and 
for small q, Re(fcg) is small. As q increases Ke{kq) 
starts to increase rapidly and Im(fcg) decreases. 
The value and magnitude of this sharp change in 
kq are determined by the proximity to ODT. Far- 
ther from the critical point (smaller x < Xc) the 
variation of kg with q are smoothed out. 

• Fig. 2: (a) A plot of the phase angle aq from Eq. 
(|l9|), as a function of the surface g-mode, for dif- 
ferent Flory parameters Nx- Circles, dotted, dia- 
mond and dashed lines correspond to Nx = 6.9, 
8.5, 9.3 and 10.2, respectively. The solid line is for 
Nx = Nxc- Far from the ODT point (high temper- 
atures, X Xc), all g-modes have phases equal to 
zero, creating a quarter-lamella region of adsorp- 
tion near the surface. In the opposite limit, i.e. 
when X ~ Xc g-modes with q < qo have aq 7r/2, 
while the g-modes with q > qo have aq — 7r/2. 

In (b) are shown the surface amplitudes |Aq + 
from Eq. ([ig|), as a function of the q-mode, for the 
same series of x values as in (a). As x ~^ Xc the 
q = qo mode goes critical first, with a diverging 
amplitude. 

• Fig. 3: A contour plot of the BCP order parame- 
ter (j)(x,y), where the surface pattern (bottom line, 
y — 0) contains only one mode: (t{x) — Uq cos{qx), 
with (jq = 1. A-rich regions are black while B-rich 
are white. In (a), q — qo while in (b) q — O.Sgo- 
The decay length ^ is smaller in (a) than in (b) 
because the surface q mode is larger. The Flory 
parameter was set to Nx = 10.4 < A^Xc- 

• Fig. 4: A contour plot of the BCP order parameter 
close to the critical point {Nx = 10.4). The surface 
pattern at y = is 17(2;) = cro+crq coa{qx), where (Tq 
is the average preference and q = |go is the modu- 
lation g-mode, with amplitudes o-q = <^q = 0.1. The 
g = surface mode has a longer range effect than 
the g > surface mode, and induces parallel lamel- 
lar arrangement farther away from the surface. At 
yet larger distances the order parameter decays to 
its bulk <j) — value. 

• Fig. 5: A sketch of a thin-film system confined be- 
tween two surfaces. The y axis is perpendicular to 
the two parallel surfaces, located at y = zkL. The 
z axis is out of the plane of the paper. 

• Fig. 6: Polymer order parameter for a system of 
homogeneous interactions = 1 at y — L — 1.5c?o 
surface and modulated interactions at the opposite 
surface a~ = cos(gx), y = —L. The surface modu- 
lation wavenumber was chosen to be g = 0.5go. As 



expected, the B polymer (shown in black) is pref- 
erentially attracted to the upper surface, while the 
bottom surface exhibits modulated adsorption pat- 
tern. This pattern propagates to the top surface. 
The Flory parameter was chosen Nx — 10.4. 

• Fig. 7: The two amplitude functions fo{y) (dashed 
hue) and fq{y) (solid hue) from Fig. 6 plotted 
against y/do- fo is negative at y = L (top uniform 
surface of Fig. 6), and fq is negative at the opposite 
modulated surface, y — ~L. Notice that the pat- 
tern at y = —L induces order at the vicinity of the 
other surface, as fq{L) ^ 0, although cr"*" = const. 

• Fig. 8: A system of modulated surface = 
cos(ga;) at y = —L = —l.bdo and of uniform = 1 
at the opposite surface, y = L = 1.5c?o, for a series 
of different values of q/qo- The effect of chang- 
ing the repeat period g is clearly seen when q/qo 
varies: 1 in (a), 2/3 in (b), 1/3 in (c). In all cases, 
Nx = 10.2. 

• Fig. 9: The crystalline- like checkerboard character 
of polymer order parameter. In (a) two patterned 
surfaces in phase with one another, but opposite 
in sign: cr"'" = —a~ = aqCos{^qox). In (b) the 
patterns are 7r/4 out of phase: = aq cos{^qox), 

= (Tg sin(igoa:). The amplitude is aq = 0.2, 
Nx = 10.4, and the top and bottom surfaces are 
located at y = ±1.25do. 

• Fig. 10: Propagation of surface order into the 
bulk, (a) is the original chemical pattern (the let- 
ters "BCP") at the y — surface, whose size is 
20 X 20 in units of do . White corresponds to A-block 
preferring regions. A sequence of contour plots for 
y — 0.5, 2do, 3.5(io: ^do and G.bdo are shown in 
(b), (c), (d), (e) and (f), respectively. The original 
pattern is gradually fading (small features, high g- 
modes first) as y is increased, until it is completely 
washed out. For y w [n+^)do with n integer, there 
is an inversion of the original pattern, as the A 
(white) and B block (black) are interchanged rela- 
tively to the original pattern. The Flory parameter 
is taken as Nx = 9.5. 

• Fig. 11: Contour plots as in Fig. 10, but in (a) the 

surface pattern is reduced to smaller size of about 
4do X 4do, while in (b) size is as in Fig. 10 but 
the temperature is higher, Nx = 9. The lamellar 
features along the letter are less prominent than 
in Fig. 10 (b). Note that in (b) the bulk ordering 
cannot tightly follow the surface pattern when the 
pattern size becomes comparable to Ido, as in (a). 

• Fig. 12: Creation of a complex three dimen- 
sional morphology by superposition of two one- 
dimensional surface patterns: a~ = cos{qx + qz) 
and cr"*" = cos(ga; — gz). Shown is the thin film BCP 
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morphology, where (a) is the surface pattern at the 
y = —L = —do surface and (c) is the pattern at 
y = L = do- A contour plot of the order parameter 
in the mid-plane, y = 0, is shown in (b), depict- 
ing A-rich and B-rich regions with square symme- 
try. The Flory parameter is taken as Nx = 9 and 
q = go/3. 

• Fig. 13: Appearance of curved lamellae as a re- 
sult of a one-dimensional surface pattern along the 
z surface axis. In (a) is surface stripe is shown in 
the X — z plane. The (white) stripe has a surface 
field of a = 0.5 inducing preferential adsorption of 
the B-polymer. The rest of the surface (denoted 
by a grey color outside the stripe) has a = and 
is indifferent to A/B adsorption, (b) is a contour 



plot in the x — y plane, depicting curved lamellae 
surrounding the "disturbance" at the middle. As 
the distance from the stripe is increased more than 
lOdo shown here, the lamellae gradually fade away. 
The Flory parameter is taken to be Nx = 10. 

• Fig. 14: Same as in Fig. 13, but with a = 0.5 in- 
side the stripe (white) , while the rest of the surface 
(black) has cr = — 0.5. In (a) the x — z surface is 
shown while in (b) the contour plots are shown in 
the X — y plane. Far from the stripe the B-polymer 
(in black) is adsorbed to the surface, and overall a 
lamellar morphology parallel to the surface is seen. 
Close to the stripe disturbance these lamellae are 
modified, distorted locally by the presence of the 
stripe. 
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